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ABSTRACT 



Serendipitously we have discovered a rare, bright submillimeter galaxy (SMG) with a flux density of 30 ± 2 mjy at A = 1.2mm, using 
MAMB02 at the IRAM 30-meter millimeter telescope. Although no optical counterpart is known for MM18423+5938, we were able 
to measure the redshift z = 3.92960 ± 0.00013 from the detection of CO lines using the IRAM Eight Mixer Receiver (EMIR). In 
addition, by collecting all available photometric data in the far-infrared and radio to constrain its spectral energy distribution, we derive 
the FIR luminosity 4.8 10 I4 /m L G and mass 6.0 l(f /m M Q for its dust, allowing for a magnification factor m caused by a probable 
gravitational lens. The corresponding star-formation rate is 8.3 10 4 /m M Q /yr. The detection of three lines of the CO rotational ladder, 
and a significant upper limit for a fourth CO line, allow us to estimate an H 2 mass of between 1.9 10"/™ M Q and 1.1 10 I2 /m M Q . 
The two lines C l( 3 Pi - 3 P ) and C l( 3 P 2 -^i) were clearly detected and yield a [Cl]/[H2] number abundance of between 1.4 10 5 and 
8.0 10~ 5 . Upper limits are presented for emission lines of HCN, HCO + , HNC, H 2 0, and of other molecules. The moderate excitation 
of the CO lines is indicative of an extended starburst, and excludes the dominance of an AGN in heating this high-z SMG. 
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1. Introduction 

Deep blank-field millimeter and submillimeter surveys of small 
fields (~ 1 deg 2 ) have revealed many dusty, starburst submil- 
limeter galaxies (SMGs) over the past decade with flux densi- 
ties of a few to about ten mjy at A — 1.2mm (e.g. Bertoldi et 
al. 2007, Greve et al. 2008), and higher at A - 850jum owing 
to dust emissivity (e.g. Smail et al 1997). Recently, the South 
Pole Telescope survey, less deep but much larger in sky coverage 
(87 deg 2 ), has found 47 brighter SMGs with flux densities be- 
tween 1 1 and 65 mjy at A = 1 .4mm (Vieira et al. 2010). Whereas 
redshifts of SMGs are crucial to study their physical properties, 
most of these dust-obscured galaxies have very faint or no opti- 
cal counterparts, making measurements of spectroscopic redshift 
extremely difficult or impossible (e.g. Smail et al. 2002). 

These dust-enshrouded star-forming galaxies are expected to 
be at high redshifts and are identified with the most massive 
galaxies assembled during an energetic early phase of galaxy 
formation. Their abundance appears to peak at z ~ 2.5 (Chapman 
et al. 2005, Wardlow et al 2010). Their star formation rate is 
prodigious at up to 10 3 M Q yr _1 , and the underlying starburst ac- 
tivity is believed to result from mergers (Blain et al. 2002). 

Lestrade et al. (2009) discovered serendipitously a rare, 
bright point source, MM18423+5938, at A - 1.2mm (30 mjy) 
by mapping 50 separate fields totalling a sky area of 0.5 deg 2 
with the MAMB02 bolometer camera (Kreysa et al at the IRAM 
30-meter millimeter telescope. Subsequently, some of us (PA 
and NS) searched but did not find local CO in the direction of 
the source suggesting not a young stellar object but an SMG 
instead, despite the Vieira et al's cumulative source count that 
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yields a chance as low as ~ 7% of finding a 30mJy SMG. 
MM18423+5938 is detected at 70 /mi but undetected at 24 /mi in 
MlPS/Spitzer images. It is in neither the 2MASS catalogue, nor 
the NVSS VLA catalogue (S \achz < 2.5mJy), and no optical 
and X-ray identifications are found in catalogues searched with 
NED (MM18423+593 is however outside the SDSS footprint). 
All these photometric data are summarized in Table Q] 

We show in this Letter that MM18423+5938 is a bright, 
high-redshift SMG. We present in Sect. 2 our IRAM/EMIR spec- 
troscopic observations of MM18423+5938 at millimeter wave- 
lengths that yielded our detections of CO and C i, and upper lim- 
its for other molecular species. In Sect. 3, we model the data to 
infer the dust and gas content of MM18423+5938 and its general 
properties. To compute distances and luminosities, we adopt the 
A-CDM concordance cosmological model, Ho = 71 km/s/Mpc, 
Q M = 0.27, and Q A = 0.73 (Hinshaw et al. 2009). 



2. Observations and data analysis 

Lestrade et al. (2009) detected MM 18423 +593 8 with a high but 
uncertain integrated flux density of 30 - 60 mjy, given that it 
was located close to the border of their MAMBO map. We reob- 
served MM18423+5938 with MAMBO at the IRAM 30-m tele- 
scope in the on-off wobbler-switching mode on 2010 January 
16 using the map coordinates (ff2ooo = 18 /, 42"'22.5 S + 0.2 s and 
£2000 = 59 o 38'30" + 2")andmeasured30 + 2mJy alA = 1.2mm. 

MM18423+5938 is located close to the border of the 
archived MlPS/Spitzer maps centered on the star GJ725AB 
(AOR 4199424). We determined the flux densities of the star 
GJ725AB and MM18423+5938 at 70 /mi and 24 /zm by means 
of aperture photometry, scaling with the stellar photospheric flux 
densities of GJ725AB predicted by the NextGen stellar atmo- 
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spheric model (Allard et al. 2001). We note a difference of 9", 
i.e. at the 3cr level, between the positions of the source in our 
MAMBO map (Lestrade et al. 2009) and the archived 70 /zm 
MIPS Spitzer map, which is not understood. 

To measure the redshift of MM18423+5938, we used the 
strategy of observation developed by Weiss et al. (2009) with the 
multi-band heterodyne receiver Eight Mixer Receiver (EMIRQ) 
at the IRAM 30-m telescope. The 3mm setup (E090) of EMIR 
provides 7.43 GHz of instantaneous, dual linear polarization 
bandwidth. The entire frequency range from 77.7 to 1 15.8 GHz 
in the 3mm band can be searched with six tunings spaced to pro- 
vide 0.5 GHz overlap. This range corresponds to < z < 0.48 
and 1 < z < 10 for the CO lines between (J=l-0) and (J=8-7). 
Observations were conducted from 2010 July 29th to August 2nd 
with precipitable waper vapor comprised between 3 and 7mm 
and with standard system temperatures of 11 OK for the E090 
setup. Data were processed with 16 units of the Wide band Line 
Multiple Autocorrelator (WILMA) providing a spectral resolu- 
tion of 2 MHz for the E090 setup. The observations were con- 
ducted in wobbler-switching mode, with a switching frequency 
of 1 Hz and an azimuthal wobbler throw of 100". Pointing and 
focus offsets were determined once every two hours and found 
to be stable. Calibration was done every 6 minutes using the 
standard hot/cold load absorber. The data were reduced with the 
CLASS software. 



http://www.iram.es/IRAMES/mainWiki/EmirforAstronomers 



Table 1. Photometry available for MM18423+5938 



Band 


S v (mJy) 


Ref 


1.4 GHz 


<2.5 


NVSS 


3 mm 


9+2.0 
-1.5 


this work 


2 mm 


9±3 


this work 


1.2 mm 


30±2 


this work 


lOOyum 


<600 


IRAS 


60yum 


<100 


IRAS 


70/um 


31±4 


Spitzer/MIPS 


24yum 


<0.6 


Spitzer/MIPS 


B V R I 


B mag > 21 


NED (DSS) 



CO(4-3) 
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Fig. 1. The three CO rotational lines detected, and the two C i 
lines. The vertical scale is T m t, in mK. 



Log Aobs (mm) 

Fig. 2. Available photometric data for MM18423+5938 (Table 
[TJ. These are superimposed upon our dust emission model (full 
red curve) based on the Milky Way dust model by Desert et al 
(1990) adapted for our source at z=3.93. It consists of three com- 
ponents ; the large grains at T</=45K, containing most of the 
mass (dotted line), the very small grains at hotter temperature 
(dot-dash), and the PAH (dashed line). 



We started to scan the whole 3mm band by integrating data 
for ~ 2 hrs for each tuning. We discovered unambiguously a line 
at 93.52 GHz after 20 minutes of integration during our third 
tuning on the second night, and continued to integrate dual po- 
larisation data for 1 .5 hr in total to reach an rms noise level of 
T m £=0.8 mK in 60 km s _1 channels (Fig[T]lefthand top panel). At 
this stage, we successively assumed that this line could be CO(l- 
0), CO(2-l), ... to calculate for each of these assumptions the 
corresponding redshift and predict the frequencies of the higher 
J transitions accessible in the 2mm band of EMIR (setup E150 
from 127 to 176 GHz). We then tuned to these frequencies with 
the E150 setup and swiftly detected a line at 140.26 GHz that 
corresponds to CO(6-5) at z = 3.92960 + 0.00013, in addition to 
the line at 93.52 GHz for CO(4-3). This identification of the CO 
transitions and determination of the redshift of MM18423+5938 
were carried out during the same night of July 30/3 1 in ~ 6 hours. 
The rest of the allocated time (15 hours) was used to search for 
CO(7-6) (detected), CO(9-8) (undetected, consequently CO(10- 
9) was not searched), and for other species, C i (two lines de- 
tected), and HCN & HNC(5-4), LiH(l-O) & HCO+(5-4), H 2 0„, 
H 2 O p , 13 CO(5-4), CS(8-7), and CS(9-8) (all undetected but in- 
teresting upper limits are discussed below, see Table [3). Two 
CO lines (5-4 and 8-7) unfortunately are in the atmospheric 2 
and H2O lines opacity domains and could not be observed. The 
spectra were of high quality, stable, and flat. Their mean levels 
measure the continuum flux densities at 2mm and 3mm (Table 
[1} owing to the excellent weather conditions (precipitable water 
vapor ~ 4 mm). 

Figure[T]displays the three CO lines detected, along with the 
two C 1 lines. Spectra were smoothed to a resolution of 30-50 
km/s. Gaussian models were fitted to the lines and the results are 
reported in Table [2] Line widths found are normal albeit small 
suggesting that we are observing a galaxy seen rather face-on. 



3. Results 

3.1. Dust emission 

The photometric data are collected in Table Q] To model these 
data in Fig. [2] we choose to use the well-established Milky Way 
dust model of Desert et al (1990) as a template. This model con- 
sists of three main components, PAH and both very small and 
large grains, and the emissivity slope is assumed to be 0=2. 
The large dust grains are dominant in mass and their tempera- 
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Table 2. Observed line parameters 



Line 


Vobs 


Tmb 


Sy 


AVpwHM 


I 


V 


L'/m/lO 10 




[GHz] 


[mK] 


[mJy] 


[kms- 1 ] 


[Jy kms~'] 


[kms-»] 


[K kms- 1 pc 2 ] 


CO(4-3) 


93.5249 


5.3 ± 0.6 


26.7 ± 3. 


175 ± 20 


4.95 ± 0.5 


0± 8 


19.7±2. 


CO(6-5) 


140.2695 


6.4 ± 1 


32.8 + 5. 


189 ± 19 


6.6 ± 0.6 


-7 ±9 


11.5±1. 


CO(7-6) 


163.6342 


4.2 ± 0.9 


21.5 ±4.6 


172 ± 27 


3.9 ± 0.5 


11 ± 11 


5.0±0.6 


CifPj-^o) 


99.8378 


1.9 + 0.6 


9.6 ± 3. 


225 ± 55 


2.3 ± 0.5 


-50 ± 24 


8.0±1.7 


Ci( 3 P 2 - 3 /' 1 ) 


164.1802 


4.2 ± 1. 


21.5 + 5. 


184 ±49 


4.2 ± 0.8 


8 ± 19 


5.3+1. 



Quoted errors are statistical errors from Gaussian fits. The systematic calibration uncertainty is 10% 
* The velocity is given relative to z=3. 929605. 



MM18433+5938 




\ 



2 4 6 8 

J(upper) 

Fig. 3. CO line flux density (points) or upper limits (arrow) mea- 
sured at the IRAM-30m, with the best-fit LVG model (dash 
line) computed with a density ne, = 10 3 cm 3 , = 45K, and 
Nco/AV = 3 10 18 cm^/kms- 1 . 



is therefore T,t=45K, n(H2)=10 3 crrT 3 , and a column density 
per velocity interval Nco/AV = 3 10 18 crrr 2 /kms _1 , as adopted 
to model the data in Figure [3] From the CO( 1 -0) line derived 
from our LVG model, we infer an H2 mass of 1.9 10 11 jm M 
with the conversion factor M(H 2 )/L'co=0.8 M /(K kms _1 pc 2 ) 
adopted for ULIRGs by Solomon et al (1997). This latter 
value yields a lower limit, while the standard conversion ratio 
for the Milky Way, 5.75 times higher, yields the upper limit 
M(H2)=1.1 10 12 /m M Q . The true mass must lie between these 
two values, and indicates large amounts of molecular gas even 
if allowing for an amplification m. Assuming a typical intrinsic 
extension of 3kpc (e.g. Tacconi et al 2006), the surface filling 
factor of the molecular component is 0.3 for m=l, and 0.03 for 
;7i=10. The Lfir/L'co of 2 10 3 L /K kms _I pc 2 is consistent 
with ratios of other luminous infrared galaxies within the scatter 
(Iono et al 2009). 



ture is estimated to be T^ = 45K, constrained by the Rayleigh- 
Jeans part of the emission. The very small grains are made of 
two temperature components, ~ 80 and ~ 130K, constrained by 
the 70^m Spitzer measurement and the robust relationship be- 
tween L C0(3 _ 2) and Lfir found by Iono et al (2009). A single 
temperature component for the small grains in the model was 
tried but the resulting Lfir is significantly inconsistent with this 
relationship. The total mass of the dust required by this model 
is Mrf I( . s; =6.0 10 9 /m M , and for a gas-to-dust mass ratio of 150, 
we infer that M ?fls =9.2 10 jm M . We suspect there is a gravi- 
tational lens along the line of sight, with an amplification factor 
m, or these values would be implausibly at least one order of 
magnitude larger than for rare hyperluminous objects. We de- 
rive the total FIR luminosity Lfir = 4.8 10 14 /m L and the star- 
formation rate SFR = 8.3 10 4 /m M /yr, by applying the relation 
of Kennicutt(1998). 

3.2. CO lines 

The CO SED of MM18423+5938 in Fig. [3] indicates moder- 
ate line excitation, peaking only at J=5; CO SEDs peak at J=6 
or 7 in local starbursts such as M82 and NGC253 (Weiss et al 
2007), peak at higher J in AGN-dominated sources, e.g. J=10 in 
APM0827, and reach a plateau for J > 8 in Mrk231 (van der 
Werf et al 2010). We ran several LVG models, to constrain the 
H2 volumic density and the kinetic temperature (cf Combes et 
al 1999). The moderate excitation implies a regime of low tem- 
perature and/or density. The gas kinetic temperature is taken to 
be equal to the dust temperature (T^=45K). When models are 
run with higher T*, they all imply n(H 2 ) lower than 10 3 ctrT 3 , 
which is not realistic for CO(7-6)-emitting clouds. Our estimate 



3.3. Atomic carbon lines 

The two lines C i( 3 / , i- 3 P ) and C 1(^2-^1) were clearly detected. 
They have comparable central velocity and line width (Tabled, 
implying that they originate from the same region in the source. 
The relation between the integrated C i( 3 Pi - 3 Po) brightness tem- 
perature and the beam averaged Ci column density with the 
usual assumption of the optically thin limit is given by 

8nkv 2 w I T it C 

/ic'Aio 3 J 

where Q(T ex ) = 1 + 3e~ r ' /7 « + 5^ T2/T " is the C 1 partition func- 
tion, and T\ — 23.6 K and T2 = 62.5 K are the energies above the 
ground state. When dealing with high-z sources, we can use the 
definition of the line luminosity (e.g. Solomon etal. 1997) and 
derive the C 1 mass via (cf Weiss et al 2003, 2005) 

Ma = 1.902 x 10- 4 e(r ex )e 23 - 6/7 » L'^^JM^]. 

The mass estimated from the higher-excitation line is expressed 
in an analogous way, and we can then deduce that 

fB^±l> =2.087 e- 38 - 9 ^, 

where line luminosities are given in Table [2] The derived 
T ex is 33.9 K. The mass of atomic carbon thus amounts to 
Mci=1.0 10 8 /m M . Given our lower and upper limits to the 
H2 mass from the CO lines, the [C i]/[H2] number abundance is 
between 1.4 10 and 8.0 10~ 5 . This is somewhat higher than the 
average [C i]/[H2] number abundance found in comparable star- 
forming objects (e.g. Barvainis et al 1997, Pety et al 2004, Weiss 
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Table 3. Line upper limits 



Line 




5 v (3cr) 


L'/10 10 




[GHz] 


[mjy] 


[K kms" 1 pc 2 ] 


v^^jyy — o ) 


1U. J'T't 


v. 1U. 




13 CO(4-3) 


89.412 


< 10. 


<8.1 


HCN& HNC(5-4) 


90. 


< 7. 


<5.5 


LiH(l-O) &HCO+(5-4) 


90. 


< 7. 


<5.5 


H 2 (1, 1,0-1,0,1) 


112.978 


< 10. 


<5.0 


H 2 0„(2,l, 1-2,0,2) 


152.554 


< 12. 


<3.3 


CS(8-7) 


79.488 


< 11. 


<11.2 


CS(9-8) 


89.420 


< 7. 


<5.6 



et al 2003, 2005, Riechers et al 2009, Danielson et al. 2010). 
In these latter estimates, although the observed Lci/Lco values 
are comparable, the [C i]/[H2] number abundances are somewhat 
dissimilar because of the various CO-to-H2 conversion factors 
adopted by these authors. Abundance lower than 1.8 10~ 5 has 
been reported (Casey et al 2010). The contribution of the atomic 
carbon to the cooling is low, Lci/Lfir = 2.5 10~ 6 , comparable to 
that of nearby star-forming galaxies (Gerin & Phillips 2000). 

3.4. Other lines 

We searched for high-density tracers, such as HCN, HNC, and 
HCO + , in particular their lowest level available, i.e. J=5-4. The 
upper limits found (Table [3]) confirm that, on average, the H2 
density is not high, as found by our LVG models of the CO line 
excitation. The HCN luminosity is higher than one third of the 
CO luminosity in local AGN-dominated objects (Imanishi et al 
2004), and we note that our 3<x upper limit (LJ 1CN(5 _ 4) /Lp 0(5 4) ~ 
0.28) is close to this limit. However, observation of HCN(l-O) 
and CO(1-0) are needed to conclude. 

We also searched for H2O emission, using the first ortho and 
para lines in their ground states, i.e. H20(/K: o /s: r = lio— >loi) and 
H20{Jfc a K c =2i 1 — >2o2) but only obtained interesting upper limits. 
Assuming that these water lines are optically thick, these upper 
limits yield a filling factor of dense clumps lower than 30% that 
of CO clouds. 

At z=2.3, a tentative detection of the H20(/a; o /s: c =2ii— >2o2) 
line was reported in IRAS F10214 (Encrenaz et al. 1993; Casoli 
et al. 1994), while at z=0.685, the fundamental transition of 
ortho-water, H20(Jk„k c = Iio— *1oi) > was detected in absorption 
towards B02 18+357° (Combes & Wiklind 1997). Water lines 
were searched for other starburst galaxies at high z (Riechers 
et al., 2006, Wagg et al.2006), and significant upper limits set. 
A search for ^0(7^/^=1 10— >loi) emission toward the z=3.91 
quasar APM 08279+5255 provided a surface-filling factor lower 
than 12% of the CO one (Wagg et al. 2006, Weiss et al. 2007). 

Finally, other interesting molecules were undetected in the 
observed bands, although with insufficient sensitivity: 13 CO(4- 
3) (limiting the 13 CO/ 12 CO emission ratio to < 1 /2), CS(8-7), 
LiH (1-0) lines, several lines of formaldehyde H2CO, and SiO, 
which is a tracer of shocks. Their lowest transitions observed are 
H 2 CO(7(2,5)-8(0,8)), H 2 CO(6(l,6) - 5(1,5)), SiO(9-8)and 
upper limits are L' < 20 10 10 K kms -1 pc 2 . 

4. Discussion and conclusion 

We have found that the brightest SMG in the North, 
MM18423+5938, is at a redshift z = 3.92960 ± 0.00013. This 
source is part of a most interesting population of similar ob- 



jects recently found by the Herschel Atlas survey (Negrello et 
al 2010 submitted). From our modelled SED, the FIR lumi- 
nosity 4.8 10 14 /m L and mass 6.0 10 9 /m M for the dust 
of MM18423+5938, and the implied star-formation rate of 
8.3 10 4 /ot M /yr, lead to suspect a gravitational lens along the 
line of sight with an amplification factor m. From our mod- 
elled CO SED, we have found that M(H 2 ) is between 1.9 and 
10.8 10 11 /m M Q depending on the conversion ratio, indicating 
large amounts of molecular gas. The CO line excitation is mod- 
erate which indicates that there is both no strong heating by a 
central AGN and a starburst that is not too extreme. The aver- 
age density of the molecular medium is low, of the order of 10 3 
cm -3 , and the gas kinetic temperature is assumed to be 45K in 
our model. The atomic carbon lines, assumed optically thin, are 
excited to T ex = 33.9K and yield M C i=1.0 10 8 /w M , which cor- 
responds to a [Ci]/[H2] number abundance of between 1.4 10~ 5 
and 8.0 10~ 5 when the limits on the H2 mass derived from the 
CO lines are used. In this high-z SMG, the C i-to-CO luminosity 
ratio is consistent with those of other high-z galaxies. 

The moderate CO line excitation found excludes a dominant 
AGN in MM18423+5938, unlike Mrk231 where CO is excited 
up to J=13 (van der Werf et al 2010). This moderate excitation 
favors an extended gas disk (typically 3kpc), rather than a com- 
pact nuclear starburst (300pc) and consequently a high CO-to-H2 
conversion ratio. This high-redshift SMG, with a star formation 
efficiency of Lfir/L'co=2400, is comparable to the lower-z sub- 
millimeter galaxies studied by Greve et al (2005). 
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